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ABSTRACT 

We investigate the star-formation properties of dynamically relaxed galaxy clusters as a function of 
cluster mass for 308 low-redshift clusters drawn from the Sloan Digital Sky Survey (SDSS) C4 cluster 
catalog. It is important to establish if cluster star-formation properties have a mass dependence 
before comparing clusters at different epochs, and here we use cluster velocity dispersion, a, as a 
measure of cluster mass. We select clusters with no significant substructure, a subset of the full 
C4 sample, so that velocity dispersion is an accurate tracer of cluster mass. We find that the total 
stellar mass, the number of star-forming galaxies, and total star-formation rate scale linearly with 
the number of member galaxies, with no residual dependence on cluster velocity dispersion. With 
the mass-dependence of cluster star-formation rates established, we compare the SDSS clusters with 
a sample of z ~ 0.75 clusters from the literature and find that on average (correcting for the mass 
growth of clusters between the two redshifts) the total Ha luminosity of the high-rcdshift clusters is 
10 times greater than that of the low-redshift clusters. This can be explained by a decline in the Ha 
luminosities of individual cluster galaxies by a factor of up to ~ 10 since z ~ 0.75. The magnitude 
of this evolution is comparable to that of field galaxies over a similar rcdshift interval, and thus the 
effect of the cluster environment on the evolution of star-forming galaxies is at most modest. Our 
results suggest that the physical mechanism driving the evolution of cluster star-formation rates is 
independent of cluster mass, at least for clusters with velocity dispersion greater than 450 km s , and 
operates over a fairly long timescale such that the star-formation rates of individual galaxies decline 
by an order of magnitude over ~7 billion years. 
Subject headings: galaxies: clusters 
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1. INTRODUCTION 

Although the inverse correlation between star- 
formation rates and galaxy density is well-established 
(e.g. Balogh et al. 1997; Lewis et al. 2002; Gomez 
et al. 2003), we do not yet understand whether this is 
due to the advanced evolution in overdense regions or to 
a direct physical effect on the star formation capability 
of galaxies in dense environments. The apparent rapid 
evolution of galaxies in dense environments (Butcher & 
Oemler 1984; Dressier et al. 1997; Poggianti et al. 1999) 
suggests that it may be possible to observe the quenching 
of star formation in dense regions by exploring clusters 
over a modest range of redshifts (0 < z < 1), though 
great care must be taken to ensure a fair comparison of 
clusters (e.g. Nakata ct al. 2005; Andrcon et al. 2006). 
This evolution may be partly related to an apparent de- 
crease in the spiral population since z ~ 1 and a cor- 
responding increase in the SO population (e.g. Dressier 
et al. 1999; Smith et al. 2005; Postman et al. 2005; 
Moran et al. 2007), yet the cause of the morphology and 
SFR evolution remains unclear. 

We have undertaken a program to directly measure 
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star-formation rates (SFRs) of cluster galaxies over this 
rcdshift range, and so far we have been able to demon- 
strate that star-formation rates of cluster galaxies de- 
pend on both rcdshift and cluster mass (Finn et al. 2004, 
2005; Poggianti et al. 2006). However, we have not 
been able to distinguish mass and evolutionary effects 
because the high-rcdshift clusters typically have lower 
masses than the low-redshift clusters. Furthermore, even 
clusters of comparable mass and redshift exhibit large 
scatter in their star-formation properties (possibly sug- 
gesting that another property of clusters, such as the in- 
tracluster medium density, may play a role; Popesso et al. 
2007; Moran et al. 2007). Thus, we need to have large 
samples of clusters with adequate mass coverage to reli- 
ably characterize cluster SFR properties. The large and 
uniform cluster catalogs selected from the Sloan Digital 
Sky Survey (SDSS; York et al. 2000) are the best avail- 
able for investigating the dependence of star-formation 
rates on cluster mass. 

Several groups have studied the environmental depen- 
dence of SDSS galaxy properties such as morphology, 
color, and star-formation rate. Often, environment is 
characterized in terms of the local density (e.g. Gomez 
et al. 2003; Blanton et al. 2005; Baldry et al. 2006), 
and it is generally found that the correlations with local 
density do not depend strongly on the larger-scale envi- 
ronment (e.g. Balogh et al. 2004; Blanton et al. 2005). 
An alternative approach is to identify bound systems and 
look for correlations with system mass, where the latter 
property is measured from galaxy dynamics (e.g. Goto 
2005; Rines et al. 2005), X-ray properties (e.g. Popesso 
et al. 2007) or via direct comparison to model predictions 
(Wcinmann ct al. 2006a). Most of these studies have fo- 
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cuscd on trying to disentangle the effects of galaxy stellar 
mass from those of environment (e.g. Kauffmann et al. 
2004; Weinmann et al. 2006a; Baldry et al. 2006). In the 
present paper our goal is to understand how the total Ha 
luminosity in clusters, integrated over all bright galaxies, 
depends on cluster velocity dispersion. This allows us to 
make a direct comparison with higher- redshift studies, 
and therefore begin to construct a star-formation his- 
tory for galaxies in clusters analogous to what has been 
done in the field (e.g. Madau et al. 1998; Hopkins 2004). 
Some attempt at this has already been made (Poggianti 
et al. 2006; Nakata et al. 2005); here we expand on these 
works, using a larger sample of local clusters with reliable 
velocity dispersions and Ha emission-line measurements. 

Specifically, we investigate star-formation as traced by 
Ha in 308 low-redshift clusters drawn from the C4 clus- 
ter catalog (Miller et al. 2005). We describe the cluster 
sample selection in $2] and the selection of Ha-emitting 
galaxies in Sj3] We quantify the mass-dependence of clus- 
ter star-formation rates in <j4j and use these results to 
quantify the environment-dependent evolution by com- 
paring with a higher-rcdshift sample of clusters in ij5l 
We summarize the main results and present conclusions 
in EJSJ We adopt a ACDM cosmology, assuming Slo = 0.3, 
Oa = 0.7, and Ho = 70 km s _1 Mpc -1 unless otherwise 
noted. 

2. CLUSTER SAMPLE & MEMBER GALAXIES 
2.1. Sample Selection 

We select our cluster sample from the C4 DR5 clus- 
ter catalog (Miller et al. 2005). We apply a minimum 
redshift cut of z = 0.05 to minimize aperture effects of 
the fiber spectra (Kewley et al. 2005) and a maximum 
redshift cut of 0.09 to minimize incompleteness of the 
galaxy catalog due to the spectroscopic magnitude limit. 
This leaves a sample of 923 clusters. 

The virialized region of a cluster scales with cluster 
mass, so we need a reliable mass estimate in order to 
properly select member galaxies for each cluster. There 
are many ways to estimate cluster mass. For example, 
Weinmann et al. (2006b) estimate halo mass using halo 
occupation models, and Rines & Diafcrio (2006) use the 
distribution of galaxies in Av — Ar space to derive clus- 
ter mass profiles. However, the simplest method, and 
the one most easily applied to higher redshift samples, 
is to use the cluster velocity dispersion. Velocity disper- 
sion can be an accurate tracer of cluster mass for relaxed 
systems that are not dominated by substructure, and for 
which there are a sufficiently large number of redshifts 
available (e.g., Miller et al. 2005). Furthermore, velocity 
dispersion can be directly related to other cluster prop- 
erties such as virial radius using the virial theorem and 
cluster size in the redshift dimension assuming an under- 
lying velocity distribution. 

We will therefore use velocity dispersion as a measure 
of cluster mass. However, for this approach to be valid a 
minimum requirement is that the clusters be in approxi- 
mate virial equilibrium, with a well-defined measurement 
of velocity dispersion. To satisfy this condition, we ex- 
clude clusters which do not appear to be dynamically 
relaxed. This selection may reduce the scatter in cluster 
properties at fixed mass compared with an unbiased sam- 
ple. It may also bias the average star-formation rates, if 
cluster mergers directly affect star formation in member 



galaxies. However the advantage is that we can more 
directly identify any correlations with system mass, in- 
dependently of dynamical effects. Clusters with no sig- 
nificant substructure should have velocity distributions 
that are nearly gaussian, and we identify unrelaxed clus- 
ters by quantifying how much their velocity distribution 
deviates from this. For each cluster, we fit a gaussian 
with a width equal to the cluster velocity dispersion from 
Miller et al. (2005). We normalize the gaussian so that 
the area within Av < 2a matches the observed number 
of cluster members in that velocity range. To illustrate, 
we show the velocity distributions and gaussian fits for 
the first twenty clusters in Figure [TJ 

We use the xi statistic to quantify the deviation of the 
velocity distribution from gaussian, where 
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Nfcn is the number of bins, N b s is the number of galax- 
ies observed in each bin, N gauss is the number of galaxies 
expected in a gaussian distribution, and o-N obg = \J Nobs 
is the error associated with the expected counts. When 
N b s is zero, we set <JN obs to 1. We use a bin size of 0.3 
times the velocity dispersion. After visual inspection of 
the velocity histograms, we find that xi < 2-0 excludes 
clusters with significant substructure within the cluster 
itself or in the nearby environment . If a cluster does not 
make the xi cu tj we try recentering the velocity distri- 
bution and then recalculate xi ■ If xt °f the recentered 
distribution meets the criteria, we keep the cluster and 
the new central velocity. 

The xt criteria is not uniform with velocity disper- 
sion; xi is systematically lower for clusters with lower 
numbers of galaxies. Therefore, we impose a minimum 
velocity dispersion of 450 km s _1 because the xi cut be- 
comes unreliable below this point, leaving a sample of 
452 clusters. In addition, we note that at all values of 
velocity dispersion this cut will preferentially select clus- 
ters with fewer galaxies, and we must keep this selection 
effect in mind when interpreting our results. 

We apply the xi < 2.0 criteria to the remaining clus- 
ters with a > 450 km s _1 , which leaves 308 clusters. We 
show velocity dispersion versus redshift for the final sam- 
ple of clusters in Figure [2l There is no trend in velocity 
dispersion versus redshift for the final sample. 



2.2. Cluster Membership 

The size of the virialized region of a cluster scales with 
cluster mass. Therefore, our radial cut for selecting clus- 
ter members must scale with cluster mass to sample anal- 
ogous regions of each cluster. We characterize the radial 
extent of each cluster in terms of R200 > which is the radius 
inside which the enclosed density is 200 times the critical 
density and approximates the virial radius of the cluster. 
Using the redshift dependence of the critical density and 
the virial mass to relate line-of-sight velocity dispersion, 
a x , to cluster mass, we express R200 a s: 



i? 2 



800 7T G p c (z) 



(2) 



Mass and Rcdshift Dependence of SF in Clusters 




-3-2-1 1 2 3 -3-2-1 1 2 3 -3-2-1 1 2 3 -3-2- 




24 




1 - 






1.6 - 


16 
8 
n 







10 12 3 





ii- 

3.4 - 


La 





-3-2-1 12 3 -3-2-1 12 3 -3-2-1 12 3 -3-2-1 12 3 

Av/ct 



Fig. 1. — Velocity histograms for first 20 clusters. The gray curve shows a gaussian with width equal to the cluster velocity dispersion 
and normalized so that the area within 2a matches the observed number of galaxies in the same velocity slice. The number in the upper 
right of each panel is the value of \i- Clusters with xi > 2.0 are rejected due to the presence of significant substructure within or near the 
cluster. 
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Fig. 2. — Cluster velocity dispersion versus redshift. The bottom 
panel shows a histogram of number of clusters versus redshift. The 
side panel shows the histogram of number of clusters versus velocity 
dispersion. 



where the critical density. p c (z), is given by the following 
equation: 
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Simplifying, we can write R200 as 
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In deriving Equation [2] we have implicitly assumed that 
galaxies arc orbiting isotropically in a single isothermal 
sphere potential, so a x is related to the circular velocity 
V c by a x = V c /V2. 

We use the mock galaxy redshift survey of Yang et al. 
(2004) to help define the optimal selection of cluster 
members in velocity space. We use the 300 h^ 1 Mpc 
simulation and select all halos with masses greater than 
10 14 M Q so that the minimum mass of the mock clusters 
is consistent with the mass threshold imposed on the 
SDSS clusters by our a > 450 km s" 1 cut. This results 
in 1062 virialized clusters containing a total of a half mil- 
lion galaxies. The velocity distributions of the mock clus- 
ters are defined to be gaussian. Thus the mock clusters 
should be comparable to our cluster sample because we 
exclude clusters whose velocity distributions deviate sig- 
nificantly from gaussian. We then select galaxies around 
each halo that have Mb, < — 18 and velocities within 
±3<t of the central velocity. In the top panels of Figurc[3l 
we show the fraction of observed galaxies that are mem- 
bers (solid line) and non-members (dotted line) in each 
radial bin as a function of projected distance from the 
cluster center, where members are those galaxies physi- 
cally located within the virial radius of the cluster. Yang 
et al. define the virial radius using Riso! their definition 
of R v ir is thus 5% larger than our value of R2001 but this 
slightly larger radial cut does not significantly affect the 
results. The left panel of Figure [3] shows completeness 



and contamination for all galaxy types, and the number 
of contaminating galaxies exceeds the number of member 
galaxies at a radius of 0.9 times the virial radius. The 
right panel shows completeness and contamination for 
late-type galaxies only, with the number of contaminat- 
ing latc-type galaxies exceeding the number of member 
late-type galaxies at a projected radius of 0.7 times the 
virial radius. In the bottom panel of Figure [3l we show 
the fraction of members and non-members per bin as 
a function of velocity cut, including all galaxies with a 
projected separation less than R V i r . For a gaussian dis- 
tribution of velocities, a 2a membership cut within the 
virial radius provides fairly high completeness (>95%), 
while increasing the velocity cut to 3er adds more non- 
members than members. When working within R2001 wc 
therefore limit our analysis to galaxies within 2a. Note 
that the results are not sensitive to the magnitude cut; 
wc find no significant difference for magnitude cuts rang- 
ing from M B j < -16 to M Bj < -19. Thus, while a 
magnitude cut of Mb, = —19 is more comparable to 
the r-band magnitude cut we impose (as discussed be- 
low), the Mb, = — 18 allows us to increase signal-to- noise 
without affecting the results. 
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Fig. 3. — The fraction of members (solid curve) and non-members 
(dashed line) per bin as a function of projected cluster-centric ra- 
dius, R p , normalized by the virial radius (top panels) and veloc- 
ity cut, At;, normalized by the cluster velocity dispersion (bottom 
panels) in the Yang et al. (2004) model. The left panels show all 
galaxies and the right panels show late-type galaxies only. The 
errorbars show the 1 — a Poisson noise associated with the galaxy 
number counts. 



Finally, we apply an absolute magnitude cut of Mr < 
—20.68 to the cluster galaxies. This corresponds to the 
spectroscopic completeness limit of r = 17 '.77 at z = 0.09, 
the maximum redshift of our cluster sample. The abso- 
lute magnitude cut takes into account the maximum k- 
correction possible (Blanton et al. 2003a). Note that this 
is only about 0.5 magnitudes fainter than the character- 
istic luminosity M* (Blanton et al. 2003a). Our analysis 
is therefore restricted to relatively bright galaxies, and 
it should be kept in mind that star formation activity 
can be substantially higher for less massive galaxies (e.g. 
Baldry et al. 2006). 

2.3. Control Sample 
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We construct a control (field) sample to compare with 
our cluster measurements and to check for systematic ef- 
fects in our analysis. For each cluster, including those 
with substructure, we center on a random galaxy within 
the DR5 redshift survey area, regardless of its magni- 
tude, and we avoid centering near clusters. We then se- 
lect surrounding galaxies with magnitudes brighter than 
Mr = -20.68 that lie within i? 2 oo and Av < 2a. We se- 
lect five such randomly positioned samples for each clus- 
ter, resulting in a control sample of 2260 independent 
pointings. We generate control samples for clusters that 
meet the redshift cuts and lower velocity dispersion limit 
cut. However, we do not limit the control sample based 
on the substructure of its parent cluster. As a result, the 
control sample spans a slightly different range of velocity 
dispersion than the cluster sample. 

2.4. Correction for Incomplete Fiber Sampling 

Not all cluster members are targeted by SDSS spec- 
troscopy, so we must correct the cluster and control sam- 
ples for incomplete spectroscopic sampling when calcu- 
lating integrated cluster properties such as the number of 
galaxies, total stellar mass, and total star-formation rate. 
To estimate the completeness of the spectroscopic sam- 
pling, wc count the number of galaxies in the photomet- 
ric catalog with an r-band magnitude greater than that 
corresponding to M r = —20.68 at the cluster redshift. 
We compare this to the number of galaxies with spec- 
tra in the same region and magnitude slice, and multiply 
the integrated cluster properties by the ratio of photo- 
metric to spectroscopic sources. For the cluster sample, 
the average completeness is 0.86 ± 0.15, 0.86 ± 0.11, and 
0.87±0.10 for radial cuts at 0.5, 1.0, and 2.0xR 200 - The 
median completeness is 0.89. For the control sample, 
the average completeness within 0.5, 1, and 2xR2oo is 
0.92 ± 0.20,0.89 ± 0.16, and 0.88 ± 0.12. The complete- 
ness estimates for the control fields are in good agree- 
ment with the expected spectroscopic completeness of 
90% (Blanton et al. 2003b), while the completeness esti- 
mates for the clusters are slightly lower due to the higher 
density of sources. 

3. STAR-FORMING GALAXIES 

3.1. Definition of Star-Forming Galaxies 

We require star-forming galaxies to have Ha luminosi- 
ties greater than 0.41 x 10 40 ergs s~ : to sample the same 
fraction of the Ha luminosity function in each cluster. 
This limit corresponds to the minimum Ha flux detected 
at z = 0.09, fH Q = 2.0 x 10 -16 ergs s _1 cm -2 , which we 
determine empirically. This luminosity limit corresponds 
to a SFR of 0.08 M yr^ 1 if we assume a typical 
extinction value of 1 magnitude at Ha (Kennicutt et al. 
1994). 

We also require star-forming galaxies to have a rest- 
frame Ha equivalent width greater than 4A, which al- 
lows us to minimize contamination from uncertainties on 
the line width measurements and stellar absorption of 
inactive galaxies (e.g. Balogh et al. 2004). However, 
imposing this cut necessarily means we will underesti- 
mate the total amount of star formation in an Ha flux- 
limited sample. To test this, we consider the subsample 
of galaxies with z < 0.075, where we can reliably push 
to lower EW limits. We find that at most 5% of the to- 
tal star formation in our flux-limited sample occurs at 



EW< 4A. As this small contribution is somewhat un- 
certain (since measurement uncertainties and stellar ab- 
sorption still play a role) we do not correct for it. 

We exclude AGN from the sample of star-forming 
galaxies using the AGN classification provided by Miller 
et al. (2003), which have been updated to include the 
DR5 galaxies. We find 22 ± 13% of galaxies by num- 
ber with significant AGN emission within Av < ±2a 
and lxR 2 oo, and they contribute an average(median) of 
21 ± 13%(29%) of the total Ha flux. This result is con- 
sistent with the AGN fraction measured by Miller et al. 
(2003) and Kauffmann et al. (2003a). 

3.2. Star- Formation Rates & Stellar Mass 

To calculate SFR from observed Ha flux, we first apply 
an aperture correction which is computed by comparing 
the flux within the i?-band petrosian radius with the flux 
in the fiber radius (Balogh et al. 2004). We then convert 
Ha flux to luminosity by multiplying by 47rd 2 , where cLl 
is the luminosity distance. We use the Kennicutt star- 
formation relation (Kennicutt et al. 1994) to convert Ha 
luminosity to SFR, where 

1 erg s" 1 = 7.9 x 10~ 42 M Q yr _1 . (5) 

We correct for extinction using values of A 2 from 
Kauffmann et al. (2003b), assuming the extinction law 
is of the form t x oc A~ 7 and thus A r = A z + 0.23. The 
Kauffmann et al. sample has been updated to include 
DR4 galaxies but not DR5. Therefore, 23% of the galax- 
ies do not have an extinction value from these catalogs. 
To estimate the extinction for these galaxies, we use the 
DR4 galaxies to define a relationship between A r and M r . 
The median value of A r for galaxies brighter than our 
magnitude cut (M r < —20.68) remains fairly constant at 
a value of 0.76 magnitudes. Therefore, we assume 0.76 
magnitudes of extinction at Ha for the galaxies not in- 
cluded in the Kauffmann et al. (2003b) catalog. The 
standard deviation in A r about the mean is 0.30 magni- 
tudes, and changing the extinction value by one standard 
deviation alters the inferred star-formation rates by 30%. 
Thus, the error associated with the extinction correction 
is small compared to other sources of error and does not 
significantly affect our results. 

Stellar mass is the integral of a galaxy's star-formation 
history and enables us to study past star-formation effi- 
ciency as a function of cluster mass. Again, we use stellar 
mass estimates from Kauffmann et al. (2003b). To es- 
timate the stellar mass for galaxies not included in the 
Kauffmann et al. catalog, we use the DR4 galaxies to 
define a relation between M r and stellar mass. We show 
the result in Figure [U where the blue line shows the rela- 
tionship we use for galaxies with M r > —21.5 (blue line) 
and the cyan line shows the slightly flatter relationship 
for galaxies with M r < —21.5. 

4. MASS DEPENDENCE OF CLUSTER SF PROPERTIES 

We first investigate how the number of cluster mem- 
bers brighter than M r = —20.68 depends on cluster ve- 
locity dispersion. In Figure [5l we show the number of 
galaxies within R200 and with Av < 2a for the clus- 
ter (black circles) and control fields (open squares). The 
small black points show the values for the individual clus- 
ters. The number of galaxies is strongly correlated with 
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Fig. 4. — Distribution of stellar mass versus M r for the DR4 
galaxies. Blue and cyan line show average relation used to estimate 
stellar mass from M r for galaxies not included in the Kauffmann 
ct al. (2003b) catalogs. The dashed green vertical line shows the 
magnitude cut used when selecting cluster members. 



velocity dispersion for both samples. Note that the num- 
ber of galaxies in the control sample increases as approx- 
imately cr 3 as expected since the cylindrical selection vol- 
ume scales in this way. The cluster fields are overdense 
relative to the field, by a factor of 6 — 23 that appears to 
depend on velocity dispersion. Note that the clusters are 
not corrected for residual field contamination; this figure 
demonstrates that we are overestimating the number of 
cluster galaxies by no more than about 10 per cent due 
to this effect. The filled squares show the control data 
rescaled by a factor of 7.8, to allow easy comparison of 
the slope of the cluster and control samples. 



• Cluster 
□ Control 




500 1000 



a (km s 

Fig. 5. — The total number of galaxies as a function of cluster 
velocity dispersion for all galaxies with velocities within ±2<x and 
projected cluster-centric distance less than 1 X i?200- The small 
black points show the values for individual clusters, and the larger 
symbols show the average in equally populated bins for the cluster 
(black circles) and control (open squares) samples. Errorbars show 
the error in the mean in each bin. The filled squares show the 
control data, scaled so that the average counts in the control fields 
equal the average counts in the cluster fields. The solid line shows 
the best power-law fit for the cluster data, and the dashed line 
shows the best power-law fit for the control data. 



It is evident that the number of galaxies in the clus- 
ter fields follows a flatter trend versus velocity dispersion 
than the number of galaxies in the control fields. How- 
ever, we note that there are several important biases here 
that must be accounted for. First, the uncertainties on 
both values are correlated, such that an overestimated 
a leads to a larger volume and, hence, an overestimated 
iVg a i. For the field sample, the uniform galaxy density 
means one expects 7V ga i oc cr 3 , so dN %a \ oc cr 2 da. How- 
ever, for the cluster population, the galaxy overdensity 
in a given cluster is expected to decline as approximately 
R~ 2 , so the total number of galaxies within i?2oo (which 
includes the cluster and surrounding field population) 
should scale like oc a a , with a < 3. At the extreme, 
where the background population is negligible and most 
of the cluster population is well within i?200j -^ ga i is in- 
dependent of uncertainties in a. For these reasons, at 
a fixed value of cr, observational uncertainties alone will 
scatter the cluster measurments of N ga \ along a slope that 
is shallower than that of the field measurements. We do 
not attempt to quantify this however, as it is sensitive 
to the radial distribution of galaxies within clusters, and 
the relative contribution of foreground and background 
galaxies. Finally, cluster selection biases will play a role. 
The imposed limit on cr results in a bias at low cr toward 
clusters where a is overestimated. Furthermore, our \ 
cut on velocity distribution tends to favor systems with 
fewer members, which will result in a bias at high a to- 
ward fewer N ga i . For all these reasons we do not attribute 
much significance to the apparently flatter slope of the 
clusters. Hereafter, when attempting to quantify cluster 
star-formation properties as a function of cr, we will nor- 
malize by the total number of member galaxies so that 
we reduce the cr-depcndcncc of the dependent variable. 
Similar results would be obtained if we normalized by 
the total stellar mass (we will show below that the two 
are strongly correlated, as expected); however for com- 
parison with high-rcdshift studies it is preferable to use 
a simple measurement like N ga i if possible. 

4.1. Distribution of Ha Luminosities 

We show the distribution of Ha luminosities as a func- 
tion of cluster velocity dispersion in Figure [H We split 
the cluster sample into three groups based on mass, and 
dividing at a =575 km s _1 and cr =700 km s" 1 yields 
an approximately equal number of galaxies per group 
(~ 800). The points (stars, open circles, open triangles) 
show the number of galaxies per luminosity bin divided 
by the total number of galaxies in that mass division 
(~ 800 galaxies), and the errorbars show the 1 — cr Pois- 
son noise. The three sets of points arc consistent within 
the errors, showing that the distribution of Ha luminosi- 
ties is not a strong function of cluster velocity dispersion. 

In Figure we compare the distribution of Ha lu- 
minosities for the cluster (stars) and control (filled cir- 
cles) samples. The open circles show the distribution for 
the control fields after being scaled to meet the average 
counts of the cluster fields. The errorbars for the clus- 
ter and scaled control fields overlap, indicating that the 
shapes of the cluster and control Ha luminosity distribu- 
tions are similar. This implies that the primary difference 
between the cluster and control sample is in the fraction 
of star-forming galaxies, with the star-forming fraction of 
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Fig. 6. — Number of galaxies per bin divided by total number of 
galaxies versus L(Ha). Stars, circles, and triangles show galaxies 
from clusters with u < 575, 575 < a < 700, and a > 700 km s _1 , 
respectively. Errorbars reflect 1 — a Poisson noise. 



the control fields a factor of ~ 1.6 higher than the cluster 
fields. This is consistent with previous results based on 
smaller samples (Carter et al. 2001; Balogh et al. 2004; 
Rines et al. 2005), which have shown that the distri- 
bution of Ha luminosities and EWs of the star-forming 
population is independent of environment . 
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Fig. 7. — Number of galaxies per bin divided by total number 
of galaxies versus L(Ho) for the cluster (star) and control (filled 
circles) samples. The open circles show the distribution for the 
control fields after being scaled to meet the average counts of the 
cluster fields. Errorbars reflect poisson noise. The 



4.2. Integrated Cluster Properties 

We now look for a mass dependence in the integrated 
properties of the star-forming cluster galaxy population. 
In Figure [3 we show the total stellar mass (panel 1), 
the number of star- forming galaxies (panel 2), and the 
total SFR (panel 3) for all galaxies within R200 and with 
Av < 2a versus cluster velocity dispersion, where all the 
dependent variables have been normalized by the num- 
ber of galaxies in each cluster to reduce our sensitivity to 



correlated errors. The left and right columns of Figure 
[S] show the cluster and control data, respectively. The 
small black points show the values for individual fields, 
and the filled circles show the average in equally popu- 
lated bins. The errorbars, often smaller than the filled 
circles, show the error in the mean. The solid line shows 
the best power-law fit in each panel, and the dashed hor- 
izontal line shows the average value to show what one 
would expect if there is no dependence on cluster mass. 
None of the power-law fits deviates significantly from a 
flat relationship, indicating all three measures of star- 
formation efficiency are independent of cluster mass, at 
least for clusters with a > 450 km s _1 . An equivalent 
way to express these results is that the total stellar mass, 
the number of star-forming galaxies, and the total SFR 
scale linearly with the number of galaxies, independent 
of cluster velocity dispersion. 

Our results are generally consistent with those of other 
studies that have measured star-formation efficiency as a 
function of cluster mass using smaller samples of SDSS 
clusters. Poggianti et al. (2006) use a smaller subset of 
the C4 clusters and find that the fraction of star-forming 
galaxies decreases as velocity dispersion increases from 
~300 to ~500 km s _1 and then levels off. We do not 
probe to low enough velocity dispersions to compare re- 
sults at a = 300 km s _1 , but our results arc consistent 
with a flat relationship for a > 450 km s . Note that 
this result is unchanged if we keep AGN in our sample 
(e.g. Poggianti et al. 2006). These results are fully con- 
sistent with the lack of mass dependence observed from 
smaller samples of both optically-selected (Goto 2005) 
and X-ray selected (Popesso et al. 2007) clusters. 

5. EVOLUTION OF CLUSTER STAR FORMATION 

We now use our SDSS cluster sample as a low-redshift 
baseline and compare with star-formation properties of 
higher-redshift clusters from the literature. We will fo- 
cus on those few studies that have comparably complete 
observations of Ha emission. 

There are a growing number of Ha surveys of galaxy 
clusters at z > 0.1, including both spectroscopic (Couch 
et al. 2001; Balogh et al. 2002) and narrow-band imag- 
ing (Balogh & Morris 2000; Finn et al. 2004; Kodama 
et al. 2004; Finn et al. 2005) surveys, which are suit- 
able for our purposes. For this comparison we use only 
the narrow-band Ha imaging data for three z ~ 0.75 
clusters from Finn et al. (2005) because they provide 
the longest redshift baseline and the narrow-band ob- 
servations for all three clusters were made in a uniform 
way using the same telescope and near-infrared camera. 
In addition, the three clusters are from the ESO Dis- 
tant Cluster Survey (EDisCS; White et al. 2005), and so 
they have readily available derived-data products (such 
as k-corrected absolute magnitudes and estimates of the 
total number of member galaxies; Pello et al. 2007; Rud- 
nick et al. 2007) that enable a more precise compari- 
son with the SDSS clusters. Furthermore, spectroscopic 
(Halliday et al. 2004; Milvang- Jensen et al. 2008), weak 
lensing (Clowe et al. 2006), and X-ray studies (John- 
son et al. 2006) show that the velocity dispersions pro- 
vide a reliable estimate of cluster mass for these clusters. 
The EDisCS clusters CL 1040-1155, CL 1054-1245, and 
CL 1216-1201 are at redshifts of 0.704, 0.748 and 0.794, 
respectively. Note that we do not have the data nec- 
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Fig. 8. — (1) Total stellar mass, (2) number of star-forming galaxies, and (3) the total SFR per total number of galaxies versus cluster 
velocity dispersion for all galaxies with velocities within ±2<r and projected cluster-centric distance less than 1 X i?200- The left and right 
panels show the cluster and control data, respectively. In each panel, the small black points show the values for individual clusters or 
control fields, and the larger symbols show the average in equally populated bins. Errorbars showing the error in the mean in each bin 
are smaller than the symbols. The solid line shows the best power-law fit for the data, and the dashed line shows the average value of the 
y-axis variable. 
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essary to identify AGN in the z ~ 0.75 clusters. As a 
result, we include SDSS galaxies with AGN emission as 
well, and hereafter we discuss Ha emission in terms of 
L(Ha) rather than SFR. 

5.1. Integrated Ha Luminosity 

The integral of the Ha luminosity function gives the to- 
tal Ha luminosity per cluster or integrated SFR if AGN 
are excluded. This quantity is analogous to the volume- 
averaged SFR that is commonly used to quantify the 
evolution of field galaxies (e.g., Madau et al. 1998; Hip- 
pelein et al. 2003; Hopkins 2004), and so we use the total 
Ha luminosity as a means of quantifying the evolution 
between the SDSS and z ~ 0.75 clusters. 

When calculating the integrated Ha luminosity, we 
apply several selection criteria to ensure that we are 
comparing complete and analogous galaxy samples in 
the higher and low-redshift clusters. First, we apply 
the same radial and velocity cuts to all clusters, includ- 
ing all galaxies within 0.5xR2oo and Av < 3a. The 
radial cut matches the areal coverage of the EDisCS 
clusters, and the velocity cut matches the filter width 
for CL 1216-1201. The filters for CL 1040-1155 and 
CL 1054—1245 correspond to a velocity width closer to 
6<t. The correction for the extra field contamination in 
these filters, which we detail in the Appendix, amounts 
to scaling the integrated L(Ha) of CL 1040—1155 and 
CL 1054-1245 by a factor of 0.81. Second, we apply 
the same rest-frame equivalent width cut of 10A to all 
star- forming galaxies, which is the minimum Ha equiv- 
alent width detected in the z ~ 0.75 narrow-band imag- 
ing surveys. Third, we include only those galaxies with 
M r < -20.68, where the values of M r for the EDisCS 
clusters are taken from Pello et al. (2007). Fourth, we 
require the SDSS galaxies to have Ha luminosities that 
lie above the flux limit of the z ~ 0.75 Ha imaging, which 
we discuss in more detail in £ 15.21 Finally, we correct the 
integrated L(Ha) of the SDSS clusters by the spectro- 
scopic completeness as described in £12.41 

The cluster version of the Madau plot (Madau et al. 
1998) needs to take into account that the integrated Ha 
luminosity scales with cluster mass. Therefore, instead 
of showing integrated Ha luminosity versus redshift, we 
show EL (Ha) versus cluster velocity dispersion in Fig- 
ure [9] and measure an offset between the low and high- 
rcdshift cluster samples. The open circles show the av- 
erage EL(Ha) in equally-populated bins versus velocity 
dispersion for the SDSS clusters, and the errorbars show 
the error in the mean. The filled points show EL (Ha) for 
each individual high-rcdshift cluster, and the errorbars 
are the uncertainty of the sum. The bold solid lines in 
Figure [9] show a slope of 2.9, which is the best-fit slope 
derived from the SDSS clusters. We assume the same 
slope for the z ~ 0.75 clusters and fit the zeropoint by 
least-squares minimization. The thin dotted lines show 
the standard deviation associated with zeropoint for both 
the SDSS and EDisCS power-law as determined by boot- 
strap resampling. The offset between the best-fit lines for 
the SDSS and z ~ 0.75 clusters corresponds to a factor of 
~26±16 increase in EL(Ha) when measured at a fixed 
velocity dispersion, with the uncertainty dominated by 
the uncertainty in the high- redshift zeropoint. Clusters 
evolve, however, and a typical cluster has roughly dou- 
bled its mass since z ~ 0.75. The gray solid lines in Fig- 



ure [9] link the z = 0.07 clusters with their progenitors at 
Z = 0.75, where we have used the semi-analytic models of 
Wechsler et al. (2002) to estimate the redshift-evolution 
of cluster mass. Taking the mass evolution into account, 
the inferred evolution for a given cluster is a median fac- 
tor of ~ 10 ± 6, with the amount of evolution increasing 
as cluster mass decreases. 

! ♦ CL1 040-1 155 ' ' ' ' 
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Fig. 9. — Integrated Ha luminosity versus cluster velocity dis- 
persion for C4 clusters (open circles) and z ~ 0.75 EDisCS clusters 
(filled symbols). The solid black lines show a slope of 2.9, which 
is the best-fit slope derived from the SDSS clusters. The dotted 
lines show the standard deviation associated with the zeropoint of 
powcrlaw fits. The solid gray lines link the z = 0.07 clusters with 
their progenitors at z = 0.75. 

To measure the average evolution in the total Ha lu- 
minosity between the SDSS and highcr-redshift clusters, 
we assume that the slope of the 'EL(Ha) — a relation is 
the same for the SDSS and higher-redshift clusters. The 
change in zero point then indicates the average evolution. 
However, the slope of the relation could also vary with 
rcdshift. For example, Poggianti ct al. (2006) present a 
detailed analysis of the star-forming population for the 
complete EDisCS sample of 17 0.4 < z < 0.8 clusters, 
and their results suggest that the trend of total star- 
formation rate versus cluster mass is flatter at higher 
rcdshift. Poggianti et al. use the [O II] emission line as 
a star-formation indicator, and here we limit our com- 
parison to the 3 EDisCS clusters that have Ha luminosi- 
ties. As a result of the small sample size, we can not 
independently measure the dependence of star formation 
on cluster mass for the high-redshift sample. However, 
because the SDSS and 3 EDisCS clusters presented here 
span the same range of velocity dispersions, the zeropoint 
calculation will yield a similar number even if the slope of 
the relation changes with redshift. Thus our measure of 
the average evolution in the total Ha luminosity between 
z = 0.07 and z ~ 0.75 should be insensitive to changes 
in the slope of the EL (if a) — a relation, although the 
interpretation of this result could be different. 

As with the analysis in £14.21 we t r Y to reduce the de- 
pendence of the dependent variable, the integrated Ha 
luminosity, on a by normalizing by the total number of 
member galaxies. For the EDisCS clusters, we use the 
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estimate of the total number of galaxies that is derived 
from the luminosity functions of Rudnick et al. (2007); 
these estimates have associated uncertainties of 30-50%. 
The results, shown in Figure [TO] indicate that the aver- 
age Ha luminosity per galaxy has increased by a factor 
of 13 ± 7 between z ~ 0.07 and z ~ 0.75. We show 
the fraction of star-forming galaxies versus redshift in 
Figure [11] The fraction of star-forming galaxies (where 
star-forming galaxies are those that meet the selection 
criteria outlined above) is a factor of 6 ± 3 higher in the 
EDisCS clusters. Poggianti et al. (2006) find that the 
fraction of star-forming galaxies varies with cluster mass 
at higher redshift. Again, we are comparing with only 
3 highcr-redshift clusters, and so we do not attempt to 
quantify a mass-dependent evolution in the star-forming 
fraction. 
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Fig. 10. — Integrated Ha luminosity divided by number of mem- 
ber galaxies versus cluster redshift for C4 clusters (open circle) and 
z ~ 0.75 EDisCS clusters (filled symbols). The point for the C4 
clusters shows the average for the 308 clusters, and the errobars 
show the standard deviation. 



5.2. Distribution of Ha Luminosities 

We now try to understand the global evolution of clus- 
ter star-formation properties in terms of the evolution of 
individual star-forming galaxies. We show the low and 
high-redshift star-forming galaxies in Figure 1121 where 
we plot L(Ha) versus M r in panel (la) for all galax- 
ies within 0.5xR2oo and Av < 3er. The black points 
show the galaxies in the 308 SDSS clusters and the filled 
stars show the galaxies in the three z ~ 0.75 clusters 
(we do not correct for the additional field contamina- 
tion in CL 1040-1155 and CL 1054-1245). The dotted 
line in Figure [12] shows the approximate flux limit of the 
z ~ 0.75 imaging, and the dashed vertical line shows 
the magnitude corresponding to the SDSS spectroscopic 
completeness limit. The solid gray line shows the approx- 
imate effect of imposing a minimum Ha EW cut of 4A 
for the SDSS galaxies. The galaxies that meet the high- 
redshift flux limit and the SDSS spectroscopic magnitude 
cut are shown in panel (lb) of Figure [T2l 

To compare the distribution of Ha luminosities, we 
consider only the star-forming galaxies that lie above the 
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Fig. 11. — Fraction of star-forming galaxies versus cluster redshift 
for C4 clusters (open circle) and z ~ 0.75 EDisCS clusters (filled 
symbols). 



high-redshift flux limit and are brighter than the SDSS 
spectroscopic limit (shown in panel (lb), Figure [T^]) so 
that we are comparing relatively complete populations in 
both cluster samples. In panel (lc) we show the number 
of galaxies per bin divided by the total number of clus- 
ter members versus L(Ha) for the SDSS galaxies (open 
circles) and EDisCS galaxies (filled stars). The high- 
redshift clusters contain a higher fraction of the star- 
forming galaxies. This could be explained by pure num- 
ber evolution, if the fraction of star-forming galaxies per 
cluster increased by a factor of ~ 10, with no change in 
the distribution of Ha luminosities of star-forming galax- 
ies. However, the SDSS L(Ha) distribution peaks at a 
lower value of L(Ha) than the EDisCS distribution, indi- 
cating that there is some evolution of the Ha luminosities 
of individual galaxies. 

As an example of evolution that is due solely to an 
increase in the average star-formation rate per galaxy 
(among the active population), we simply scale the Ha 
luminosities of individual SDSS galaxies by a factor of 
10, which is the amount of evolution observed in the in- 
tegrated Ha luminosity. We show the scaled SDSS galax- 
ies in panel (2a) of Figure [12] the comparable samples of 
SDSS and EDisCS galaxies in panel (2b), and the cor- 
responding distribution of Ha luminosities in panel (2c). 
This simple scaling brings the number of star-forming 
galaxies in the SDSS clusters into agreement with the 
EDisCS distribution. Note that because we have ob- 
served only three high-redshift clusters, we do not ex- 
pect to detect high-redshift galaxies in bins where the 
expected counts fall below 0.33 galaxies per cluster. 

While this simple model does an excellent job of re- 
producing the evolution of the distribution of Ha lumi- 
nosities and is consistent with the evolution of the in- 
tegrated Ha luminosity, we must also consider that the 
broad-band magnitudes of galaxies evolve. On average, 
galaxies have faded since z ~ 0.75, and this impacts the 
inferred L(Ha) distribution because galaxies in the high- 
redshift sample may have faded below the magnitude cut 
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Fig. 12. — (la) L(H«) versus M r for SDSS galaxies from 308 clusters (black dots) and z ~ 0.75 galaxies (blue stars) from 3 clusters. 
Vertical dashed line shows spectroscopic completeness limit of SDSS. The solid line shows the EW limit of the SDSS spectra. The dotted 
line shows flux limit of z ~ 0.75 surveys, (lb) Same as panel la but showing only those galaxies that lie above the flux and magnitude 
limits, (lc) Distribution of L(Ha) for galaxies in panel lb. (2a) L(Ho) versus M r for SDSS galaxies after L(Ha) is scaled by a factor of 
10. (2b) Same as panel 2a but showing only those galaxies that lie above the flux and magnitude limits. (2c) Distribution of L(Ha) for 
galaxies in panel 2b. (3a) L(Hct) versus M r for SDSS galaxies after L(Ho) is scaled by a factor of 3 and M r is brightened by 1 magnitude. 
The red dashed line shows the magnitude of the SDSS spectroscopic completeness limit, brightened by 1 magnitude. (3b) Same as 3a but 
showing showing only those galaxies that lie above the flux and magnitude limits. (3c) Distribution of L(Ha) for galaxies in panel 3b. 



in the low-redshift sample. We account for this in a sec- 
ond model by brightening each SDSS galaxy by 1 magni- 
tude to account for fading since z ~ 0.75 (e.g. Poggianti 
& Barbaro 1997) while simultaneously scaling the Ha lu- 
minosities of SDSS galaxies by a factor of 3. The scaled 
SDSS galaxies are shown in panel (3a), and the galaxies 
that meet the EDisCS flux cut and the SDSS magnitude 
cut (which is now a magnitude brighter and shown with 
the red dashed vertical line) are shown in panel (3b). The 
resulting L(Ha) distributions (shown in panel (3c)) are 
in reasonable agreement, although the sample of EDisCS 
galaxies used to derive the z ~ 0.75 L(Ha) distribution 
is uncomfortably small. 

These models show that the evolution of the star- 
forming cluster galaxies between z ~ 0.75 and z ~ 0.07 
can be characterized by a fading of both the Ha and 
broad-band luminosities of individual galaxies. How- 
ever, the two scenarios we present in Figure [12] arc by 
no means unique, and further study is required to char- 



acterize the evolution of the Ha luminosity distribution 
in terms of the evolution of the number, magnitudes, 
and star-formation rates of the star- forming galaxies. 
Nonetheless, these preliminary results suggest that the 
Ha luminosities of individual cluster galaxies have faded 
by a factor of up to ~ 10 since z ~ 0.75. This factor 
of rsj 10 is an upper limit on the evolution, and it will 
be lower if there is evolution of galaxy magnitudes. In 
addition, galaxies that are presently on the red sequence 
may have been forming stars at z ~ 0.75 (e.g. Bell et al. 
2007), and this may further decrease the inferred evolu- 
tion of individual Ha luminosities. 

Several issues compromise our comparison of the SDSS 
and z ~ 0.75 star- forming galaxies. The dominant limi- 
tation is the small number of high-redshift clusters. We 
have shown in 21 that ESFR exhibits large scatter as a 
function of cluster mass, so a more precise measure of 
evolution requires a much larger sample of high-redshift 
clusters. A second limitation is our use of Ha as a star- 
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formation indicator. Although relatively insensitive to 
dust compared with other commonly used optical indi- 
cators such as [Oil] emission (e.g. Nakata et al. 2005), 
it is well known that a full census of star formation must 
account for energy re-radiated by dust (e.g. Due et al. 
2002; Coia et al. 2005; Bell et al. 2007; Bai et al. 
2007). If dust obscuration is a strong function of red- 
shift and environment, this could affect our conclusions. 
Furthermore, for the high-redshift clusters we are not 
able to quantify the AGN contamination through the Ha 
emission alone. AGN account for an average (median) 
of 23% (34%) of total Ha emission for galaxies within 
0.5R 200 and Aw < 3ct for the SDSS clusters and 20% 
(average) by number. Existing studies of high-redshift 
clusters show low AGN fractions (e.g., Homeicr et al. 
2005), so AGN contamination is probably not a limit- 
ing factor in the comparison. Upcoming studies of the 
high-redshift clusters at X-ray and infrared wavelengths 
will provide some insight into the AGN fraction of the 
EDisCS clusters. A final caveat that may affect our find- 
ings in 2] and $3] is that the SDSS Ha luminosities are 
derived from aperture fluxes. Koopmann et al. (2006) 
find radial gradients in Ha emission that are not reflected 
in the radial distribution of stellar light, and this com- 
promises the aperture correction we apply to the Ha flux. 
As a result, our analysis is not sensitive to processes that 
preferentially affect the perimeter of galaxies. 

5.3. Comparison with Field Evolution 

Field surveys show a factor of ~10 decrease in 
SSFR/Mpc 3 since z ~ 0.8 (e.g. Tresse & Maddox 1998; 
Glazcbrook et al. 1999; Tresse et al. 2002; Hippelein 
et al. 2003; Bell ct al. 2007). In comparison, Figure [T2l 
shows that we can reproduce the evolution between the 
star-forming galaxies in the low and high-redshift clusters 
by scaling the Ha luminosity of the SDSS galaxies by a 
factor of ~ 10. This is remarkably similar to what is seen 
in the field, and suggests that the strong evolution ob- 
served is not driven by cluster-specific physics. However, 
it is important to distinguish the average star-formation 
rate per galaxy and the average star-formation rate of 
the star-forming galaxies. Recent work by Bell et al. 
(2007), based on the COMBO-17 photometric redshift 
survey, suggests that the field evolution is due in roughly 
equal parts to an increase in the number density of star- 
forming galaxies, and in the average star-formation rate 
per galaxy among this population. Our analysis of the 
Ha luminosity distribution suggests a similar division in 
clusters, with an increase in galaxy star-formation rates 
contributing about half of the factor of ~ 10 observed 
evolution in the integrated Ha luminosity. However, we 
require a larger sample of high-redshift clusters, prefer- 
ably with complete redshifts, and a comparable sample 
of high-redshift field galaxies to confirm this. 

An issue that complicates the comparison in evolution- 
ary rates between clusters and the general field is that 
clusters have approximately doubled their mass over the 
time range we are probing. However, many of the ac- 
creted galaxies may already be in the high-redshift survey 
volume (the ±3ct velocity cut corresponds to a large vol- 
ume). A careful comparison with simulations is needed 
to understand how infalling galaxies affect our measure 
of cluster star-formation properties. 



5.4. Scenarios for Declining SFRs 

Our main result is that the rate of evolution of clusters 
does not appear to be significantly different from that of 
the field. This suggests the cluster environment has not 
had a direct effect on galaxy star-formation rates over the 
past few billion years. Either the difference we observe 
between cluster and field populations was imprinted at 
high redshift, or it is driven by physics occurring on group 
scales (or both). Here wc discuss the implications of 
our results on the various mechanisms frequently used to 
explain the decline of cluster star-formation rates. We 
remind the reader that our observations are restricted to 
relatively massive galaxies, and lower mass galaxies may 
behave quite differently. 

If dense environments arc to directly influence star- 
formation rates at low redshift (z < 1), they must pre- 
sumably do so by removing gas from galaxies, since this is 
the dominant reservoir of baryons. Typically, discussion 
of such processes is divided into three types: the rapid 
consumption of gas, through massive starbursts; the re- 
moval of cold gas bound to the disk, via ram-pressure 
stripping (Gunn & Gott 1972) or other thermal mecha- 
nisms (Nulsen 1982); and the removal of hot, more loosely 
bound gas, through similar mechanisms (Larson et al. 
1980; Kawata & Mulchaey 2007). All of these mecha- 
nisms will lead to a reduction of star formation, but they 
operate on different timescales and in different environ- 
ments. The efficacy of ram pressure stripping depends on 
the velocity of the galaxy relative the ICM and the den- 
sity of the ICM. Both factors reach their maximum val- 
ues at the cluster center, and so ram-pressure stripping 
should be most effective for galaxies plunging through 
the middle of the cluster. One also expects that ram- 
pressure stripping is more effective in higher-mass clus- 
ters because they have higher-density ICMs (e.g. Mohr 
et al. 2000; Vikhlinin et al. 2006), and their velocity dis- 
persions are higher. It is actually difficult to strip away 
all the cold gas in a disk, except in the centers of the 
most massive clusters (e.g. Quilis et al. 2000; Roediger 
& Brueggen 2007). The fact that we see no trends in 
average star-formation rate as a function of cluster mass 
suggests 7 that ram-pressure stripping of cold gas is not 
the dominant mechanism driving the evolution of star- 
forming cluster galaxies although it certainly occurs in 
some cases (Kenney et al. 2004; Koopmann et al. 2006). 
Instead, the decline in the star-formation rates of cluster 
galaxies since z ~ 1 must be primarily driven by phys- 
ical processes that do not depend on cluster mass, at 
least for clusters with velocity dispersions greater than 
450 km s _1 . A similar argument holds for the consump- 
tion of gas by starbursts, since such bursts arc most easily 
caused by galaxy mergers which should be rare in cluster 
environments (but could be common in the outskirts, or 
in groups). 

In hierarchical clustering models, clusters form in re- 
gions of high overdensity and grow through the accre- 
tion of field galaxies and groups. Most models of galaxy 
formation include a crude treatment of environmental 

7 A possible caveat is that, as noted above, Koopmann et al. 
(2006) show that star-formation is truncated at the outer edge of 
cluster spirals in Virgo, and the SDSS fiber-derived Ho fluxes are 
not able to detect processes that preferentially affect the outer radii 
of galaxies. 
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effects, where satellite galaxies of a halo of any mass 
are instantaneously stripped of their hot gas, prevent- 
ing further growth of the cold gas reservoir (e.g., White 
& Rees 1978; Kauffmann 1995; Bower et al. 2006; Cro- 
ton et al. 2006; De Lucia et al. 2006). This "starva- 
tion" mechanism (originally proposed by Larson et al. 
1980) serves to shut down star formation on ~l-2 Gyr 
timescales, starting from the instant a galaxy becomes a 
satellite in a larger halo. Undoubtedly this gas is easier 
to remove than the cold gas considered above (Kawata 
& Mulchaey 2007); however, it is becoming increasingly 
clear that the simple assumption of instant removal of 
gas on any satellite, regardless of halo mass, is too effec- 
tive at suppressing star formation (e.g. Poggianti et al. 
2006; Weinmann et al. 2006b; McCarthy et al. 2007). 
Our results also suggest that star-formation is not shut 
off immediately; rather, we see a relatively gradual fad- 
ing in L(Ha) of a factor ~ 3 — 10 over a timespan of 
~7 billion years. Starvation that operates over longer 
timescales can also explain why field galaxies evolve at 
a similar rate, where the decline in field star-formation 
rates is due to the build-up of groups. 

Another way to see this is to directly compare the 
amount of star formation in our clusters with predictions 
from the models. For this purpose, we use the publicly 
available Bower et al. (2006) models, based on the Mil- 
lennium simulation. For each parent dark matter halo 
we calculate the expected value of a and -R200 using the 
same relations used here (Eqn. [5]), and select galaxies 
based on their projected position, relative velocity (in- 
cluding Hubble flow and peculiar velocity components) 
and r-band luminosity. Using this membership criteria 
we compute the fraction of galaxies with Ha emission 
satisfying our observational criteria on flux and EW, as 
a function of velocity dispersion. The results are shown 
in Figure 1131 As in the observations, there is no signif- 
icant trend with velocity dispersion. However the aver- 
age fraction of galaxies with Ha emission is only 10%, 
significantly lower than our observed value of 20% (Fig- 
ure [5]) . This difference alone could simply be reflecting 
uncertainties in modeling Ha emission, which depends 
on the dust model, metallicity, and AGN contribution. 
But interestingly the average emission line fraction for 
the model field population is 36.5%, in very good agree- 
ment with our control sample value. Thus, the difference 
in emission-line fraction in clusters relative to the field 
is significantly larger in the models than in our data, 
and this should be robust to modeling of L(Ha) (but 
of course sensitive to environmental dust or metallicity 
effects). This lends support to our suggestion that the 
truncation of star formation in clusters is too severe in 
the current models. Note that the models assume a sup- 
pression of star formation that is not only rapid but also 
efficient in low mass systems, well below the mass range 
our cluster sample. This is because the starvation model 
as implemented has little dependence on the mass of the 
host halo, which is probably unrealistic (McCarthy et al. 
2007). Recent observations suggest that groups have sub- 
stantially higher star-formation rates than clusters (e.g. 
Wilman et al. 2005; Balogh et al. 2007), in contrast to 
these simple model predictions, and therefore these low 
mass systems are an important environment worthy of 
further study. 
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Fig. 13. — The predicted fraction of Ho-emitting galaxies ver- 
sus cluster velocity dispersion for z = clusters selected from the 
model of Bower et al. (2006). Each cluster is shown as an indi- 
vidual dot. The many clusters with no emission-line galaxies are 
represented with a mean of 0.05 and artificial scatter, for visualiza- 
tion purposes only. The galaxy magnitude limit, and Ha flux and 
EW limits, are the same as in our SDSS sample. The horizontal 
dotted line represents the average fraction among all z = galaxies 
in the model (above the same magnitude limits). The thick solid 
line represents the average emission line fraction in bins of velocity 
dispersion. 



6. SUMMARY & CONCLUSIONS 

Using 308 dynamically relaxed, low-redshift clusters 
from the C4 cluster catalog, we investigate the depen- 
dence of the total stellar mass, number of star-forming 
galaxies, and total star-formation rate on cluster veloc- 
ity dispersion. We find that the total stellar mass, the 
number of star- forming galaxies, and total star- formation 
rate scale linearly with the number of galaxies with 
no dependence on velocity dispersion. In other words, 
SL(Ha)/N ga i and Ngp/Ngai are independent of veloc- 
ity dispersion. We interpret this to mean that the clus- 
ter environment at low redshift does not differentially 
affect the star-forming properties of luminous galaxies 
(M r < —20.68) across dense environments characterized 
by a > 450 km s" 1 . 

We compare the SDSS clusters with a sample of 
z ~ 0.75 clusters from the literature and find that the 
star-forming cluster population has declined significantly 
since z ~ 0.75. Using the relationship between SLh« — cr 
defined from the local cluster sample, we measure the 
mean shift in this relationship to the high-rcdshift clus- 
ter sample. On average (correcting for the mass growth 
of clusters between the two redshifts) we find that the 
total Ha luminosity of the high-rcdshift clusters is ~ 10 
times greater than that of the low-redshift clusters. In 
more detail, we find that we are able to bring the SDSS 
Ha luminosity distribution into agreement with that of 
the high- redshift clusters by scaling the L(Ha) of individ- 
ual SDSS galaxies by up to a factor of up to ~ 10, and we 
can reduce the required evolution in Ha luminosities to 
a factor of ~ 3 if we allow for fading of galaxies between 
z ~ 0.75 and z = 0.07. Thus, using the SDSS clusters, 
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we are able to break the degeneracy between rcdshift 
and mass-evolution that limited our previous attempts 
to quantify the evolution of cluster star-formation prop- 
erties (Finn ct al. 2004, 2005). 

We compare the cluster results to the evolution seen 
in field galaxies over a similar rcdshift interval. We find 
a factor of ~10 decline in total star formation with de- 
creasing redshift for the clusters as compared to a factor 
of ~10 decrease for the field. Thus, we find no evidence 
for any differential evolution between cluster and field 
galaxies over this redshift range. This strengthens our 
conclusion that the effect of the cluster (a > 450 km s _1 ) 
environment on bright, star- forming galaxies is relatively 
modest for < z < 0.9. In addition, the mechanism 
causing the decline in galaxy star-formation rates must 
operate over a fairly long timescale such that the star- 
formation rates of individual galaxies decline by an order 
of magnitude over ^7 billion years, a decline timescale 
which is longer than galaxy evolution models commonly 
assume. 
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APPENDIX 

To estimate the additional field contamination in our high redshift surveys, due to the large velocity range sampled 
by the narrow-band filters, we use the mock redshift survey of Yang et al. (2004). We use all halos with masses 
greater than 10 14 M Q , which corresponds roughly to the mass threshold imposed by our rma > 450 km s _1 cut, 
and select all galaxies within 0.5 x R V i r and Mb, < —18. Here, Yang et al. use i?iso to define the virialized region 
of each halo; this agrees to within 5% with R2oo(the radial cut we have used throughout), and this slightly different 
radial cut does not affect the inferred contamination. We define member galaxies as those physically located within 
the virial radius and contamination as the number of non-members normalized by the number of observed galaxies. 
We calculate the contamination for late-type galaxies, which arc likely to be the ones contributing to Ha emission and 
which also have a higher contamination rate than early-type galaxies. The simulation size is 300 h' 1 Mpc, and the 
largest velocity dispersion is 2100 km/s. A ±6cr velocity cut corresponds to a volume of 252 h^ 1 Mpc, so we are not 
limited by edge effects. We find that the contamination of late-type galaxies is 25% in a ±3tr velocity cut and 31% in 
6cr cut. Therefore, we correct the integrated L(Ha) of CL 1040-1155 and CL 1054-1245 by a factor of 31/25 = 0.81 
to account for additional field contamination. The Yang et al. (2004) simulation is not an ideal comparison because 
it is at a rcdshift of zero, has B-band magnitudes, and is normalized to the 2dFGRS survey. Nonetheless, it gives a 
useable estimate of the extra field contamination for CL 1040—1155 and CL 1054—1245. 
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